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The dynamical critical exponent z is a fundamental quantity in characterizing quantum criticality, 
and it is well known that the presence of dissipation in a quantum model has significant impact 
on the value of z. Studying quantum Ising spin models using Monte Carlo methods, we estimate 
the dynamical critical exponent z and the correlation length exponent u for different forms of 
dissipation. For a two-dimensional quantum Ising model with Ohmic site dissipation, we find 2 « 2 
as for the corresponding one-dimensional case, whereas for a one-dimensional quantum Ising model 
with Ohmic bond dissipation we obtain the estimate z ~ 1. 
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I. INTRODUCTION 

Conventionally, quantum criticality can be described 
by a quantum-to-classical mapping,^] whereby a d- 
dimensional quantum model is represented by a (d+l)- 
dimensional classical model in which the extra dimension 
corresponds to imaginary time, r. It is well known since 
the work of Hert^ that this temporal dimension and the 
spatial dimensions do not necessarily appear on an equal 
footing. In the presence of dissipative terms in the ac- 
tion for instance, long-range in tera ctions are introduced 
in the imaginary time direction,!^ making the model be- 
have as if it were ((i-t-2:)-dimensional rather than {d + 1)- 
dimensional. The dynamical critical exponent z can be 
regarded as a measure of the anisotropy between the tem- 
poral dimension and the spatial dimensions, as defined by 
the scaling of the temporal correlation length, ^ C^- 
Here, £^ ^ \K — i^d " is the spatial correlation length 
upon approaching a quantum critical point K = Kc, 
with K being some arbitrary (non-thermal) coupling con- 
stant. Knowing the value of z is therefore of fundamen- 
tal importance in the study of quantum phase transi- 
tions, especially since this critical exponent determines 
the appearance of the quantum critical regime at finite 
temperatures above the quantum critical point Such 
quantum critical points with an accompanying quantum 
critical region have been suggested to be responsible for 
instance for the anomalous behavio r of the normal phase 
of high-Tc cuprate superconductors.'^ 

To illustrate the effect of dissipation on the dynamical 
critical exponent, consider first a generic (/•''-type non- 
dissipative quantum field theory. The bare inverse prop- 
agator can be obtained from the quadratic part of the 
action as + , meaning that one has isotropic scal- 
ing between the spatial dimensions and the temporal di- 
mension, i.e. z = 1. Adding local Ohmic dissipation by 
coupling each spin to a bath of harmonic oscillators,'^ the 
inverse propagator is modified to + + Assum- 
ing a phase transition to an ordered state and taking the 
limit q — > 0, w — > 0, the dissipative term \uj\ will always 
dominate over the dynamic term w^, and so, by using 
a; g^, we may naively make the prediction z = 2. Note 



that according to this argument, the dynamical critical 
exponent for a given action is independent of the spatial 
dimensionality of the system. We will refer to these scal- 
ing arguments as naive scaling, and postpone any discus- 
sion of caveats and other possible scaling choices to Sec. 

113 

If one replaces this Ohmic site dissipation with dissipa- 
tion that also couples in space and not just in time, this 
situation may change significantly. A common form of 
dissipation in the context of arrays of resistively shunted 
Josephson junctions and related models, is the Ohmic 
dissipation of gradients, i. e. of the bond variable that is 
the difference of the quantum phase between the super- 
conducting elements.!^ In Fourier space this bond dissipa- 
tion corresponds to an inverse propagator q^-|~a;^ + q^|w|. 
(See, however. Sec. IV ) Once again letting q — > 0, 



a; — > 0, we can from naive scaling expect the dissipa- 
tion to be weaker than in the onsite case since in this 
limit q^|a;| <C q^ for any positive z. A possible value is 
therefore z = 1, for which the spatial term balances the 
dynamic term and dissipation can be considered pertur- 
batively irrelevant in renormalization group sense. 

Simple arguments of the kind given above have been 
the approach most commonly used whenever a dynamical 
critical exponent is to be determined. In recent years 
there has however been progress towards calculating the 
corrections to these lowest-order estimates for z both by 
field-theoretical renormalization group methodsPHUland 
by Monte Carlo methods-'^^^'^In addition, there has also 
been considerable recent interest in dissipative systems 
exhibiting more exotic forms of quantum criticality where 
the critical exponents are varying continuously.'^SK^ 

The most notable advance from our point of view is 
however the work by Werner et al^^ justifying numeri- 
cally the naive scaling estimate for the Ising spin chain 
with site dissipation by extensive Monte Carlo simula- 
tions. More precisely, it was found that the dynamical 
critical exponent was universal and satisfied z = 2 — ry, 
with an anomalous scaling dimension 77 « 0.015. Apart 
from Ref. TJ", almost no Monte Carlo simulations have 
been performed on extended quantum dissipative models. 
(See, however, Refs . ITSl and ITOl for reviews of Monte Carlo 
simulation for dissipative systems and quantum phase 
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transitions.) The present work can therefore be regarded 
as a natural extension of the work done by Werner et ai, 
but more importantly as a first step towards more com- 
plex dissipative quantum models with bond dissipation. 
For instance, the dissipative XY model with bond dis- 
sipation is very interesting both as a model of granular 
superconductors or other systems which may be modeled 
as Josephson junction arrays.^ In particular, such a dissi- 
pative XY modeP2l and related Ashkin- Teller model^^HH] 
have been proposed to describe quantum critical fluctu- 
ations of loop-current order in cuprate superconductors. 

Finding a value of z is also of considerable interest for 
purely classical models that include strongly anisotropic 
interactions. The reason is simply that performing a 
finite-size analysis to find the critical coupling or critical 
exponents requires a choice of system sizes that reflects 
an anisotropy in the scaling of the correlation lengths. In 
other words, one ideally needs to know the relative corre- 
lation length exponent Vt/v — z a. priori for the finite-size 
analysis to be correct. 

In this work, we seek to employ Monte Carlo simula- 
tions of Ising models to answer the following questions: 
1) Can we confirm numerically that the dynamical crit- 
ical exponent is indeed independent of dimensionality? 
(neglicting the assumed small 77) 2) How will the dynam- 
ical critical exponent for Ising variables change if one re- 
places the site dissipation with dissipation that also acts 
in space? The first question will be addressed in Sec. [llj 
where we study the 2D quantum Ising model with site 
dissipation. In Sec. |III| we turn to the second question 
by studying a ID quantum Ising chain with bond dissi- 
pation in a similar manner. The results will be related 
to the naive scaling arguments for z, after which we con- 
clude in Sec. El 



II. 2D QUANTUM ISING MODEL WITH SITE 
DISSIPATION 



We first consider a quantum Ising spin model in 
two spatial dimensions coupled to a bath of harmonic 
oscillators,'^ i.e. a higher-dimensional version of the 
model considered in Ref. [131 In Fourier space, the 
quadratic part of the action for this model can be written 
as 



where a is the Ising field. The discretized real space 
representation on a L x L x LT--lattice then reads 
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We have assumed a spatially isotropic system, so that 
Kx = Ky = K. Periodic boundary conditions are im- 
plicit in the imaginary time direction and are also applied 
for the spatial directions. Note that our representation 
is equivalent to that of Ref. il3j although superficially 
appearing slightly different. 

We could, as Werner et al. , take a quantum Ising model 
in a transverse magnetic field as a starting point, and the 
field would then give rise to the quantum dynamics of the 
spins as represented by the second line in the action in 
Eq. ([2]). However, in this work we are not interested in 
the effect of a transverse field per se, and will therefore 
treat the dynamic term as a phenomenological term of 
unspecified origin. (See, however. Sec. IV I. In the fol- 
lowing, we will fix the value of the dynamic coupling of 
the Ising field to Kr = -l/21n(tanhl) sa 0.1362 and 
vary the spatial coupling K. For the (l-l-l)-dimensional 
modeP^ this choice ensures that Kc = 1 for a = 0, 
whereas in the d = 2 case it is chosen primarily for com- 
putational convenience, and to allow for direct compar- 
ison with the d = 1 case. For the Monte Carlo sim- 
ulations we have used an extension of the Wolff cluster 
algorithnPUby Luijten and Blote^ which very effectively 
treats the long-range interaction in the imaginary time 
direction. We have mainly used an implementation of the 
Mersenne Twistei^ random number generator (RNG), 
but also confirmed that other RNGs yielded consistent 
result. We also make use of Ferrenberg-SwendserP^ 
reweighting techniques which enable us to vary K con- 
tinuously after the simulations have been performed. 

We will first present the phase diagram for this model 
in the a-K plane, as shown in Fig. [l] The phase diagram 
for the (2+l)-dimensional model is very similar to that 
for its (l-l-l)-dimensional counterpart, with a disorder- 
order phase transition for increasing dissipation and/or 
spatial coupling. Along the a-axis, a temporally ordered 
state is reached at a = ac through a purely dissipative 
phase transition when K = 0, in which case the model is 
simply a collection of decoupled (0+l)-dimensional dis- 
sipative two-level systems. The long-range interaction 
in the temporal chains decays as l/jr^ — Tjp, accord- 
ingly, the phase transition is of a kind closely related to 
the Kosterlitz-Thouless transition,'^ in which the ordered 
phase consists of tightly bound kinks and antikinks. 

With the same temporal coupling values as for the 
d = 1 case we can with relative ease determine the 
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FIG. 1: Phase diagram for the 2D quantum Ising model with 
site dissipation for Kt = — 1/2 in (tanh 1). The ordered phase 
is found for large values of spatial coupling K and dissipation 
strength a. 



critical dissipation strength ac for the independent 
subsystems, see the result stated in Ref. [T31 

We have chosen a somewhat more quantitative ap- 
proach to determine the dynamical critical exponent z 
than the one given in the presentation of Werner et ai, 
so we will use the exposition of our results to detail the 
method. This method is essentially the same as the one 
applied by the authors of Refs. 1^ and for spin glasses 
in a transverse field, but as it is rather scantily described 
in the literature, we include it here for completeness. 

The basis of our approach is as follows. For systems 
with isotropic scaling, a well known method to determine 
the value of the critical coupling is to calculate the Binder 
ratio 



and use this to plot the Binder cumulant g = I — Q/3 
as a function of coupling for several {e.g., cubic, in the 
(2+l)-dimensional case) system sizes. The Binder cu- 
mulant at the critical coupling is independent of sys- 
tem size (to leading order in L), and the crossing point 
of g{K) for two different system sizes thus defines the 
(pseudo)critical point. 

However, this finite-size scaling approach breaks down 
when the system size scales anisotropically. In this case 
the scaling at criticality is given as a function with two 
independent scaling variables instead of just one, 

QiL,Lr)^g(^^,^y (4) 

and anisotropic systems according to oc are the 
appropriate choice instead of cubic systems. Hence, given 
the value of z, one should also observe data collapse as a 
function of L^./ for the Binder cumulant curves at the 
critical point. 



In order to find z self-consistently, we consider first 
the Binder cumulant as a function of Lt for given a, K 
and L. For very small L^., the system appears effectively 
two-dimensional, and consequently the increased influ- 
ence of fluctuations makes this system more disordered 
than the corresponding three-dimensional system. In the 
opposite limit of L,- — c» the system appears effectively 
one-dimensional, and with Lr ^ ^t- the system is again 
disordered. As g is a measure of the degree of order in the 
system, g — > in both the above limits, and accordingly 
g must have a maximum for some finite value L^. — L*^. 
One way of interpreting L* is as the temporal size for 
which the system appears as isotropic as it possibly can 
be (or optimally three-dimensional), the anisotropic in- 
teractions taken into account. 

The details of our procedure are as follows. First, we 
sample the Binder ratio as a function of coupling K for 
a large number of system sizes. For each value of L we 
choose at least 14 values of Lr close to the presumed 
peak position L* for the extent of the imaginary time di- 
mension. The procedure for estimating z then follows in 
three steps. For each K, curves of the Binder cumulant 
g for all L are plotted as a function of Lr, corresponding 
to the plot shown in panel (a) of Fig. |2] Second, a 4th 
order polynomial fit is made to these curves, localizing 
the points {L*,g*) defining the peaks of the functions 
giL-r) with good precision. The obtained values for the 
peak Binder cumulants for each L are then plotted as a 
function of K, as shown in panel (b) of Fig. [2] A value 
for the critical coupling Kc can be found by estimating 
the value K to which the crossing point for two subse- 
quent values of L converges for 1/L -> 0. The third step 
for finding the dynamical critical exponent is a simple 
finite size scaling analysis of the peak positions L* of the 
curves g{Lr) as shown in panel (c) of Fig. [2j assum- 
ing the relation L* — aL^, with a being a non- universal 
prefactor. Finally, one may check the self-consistency of 
the obtained values for and z by plotting the puta- 
tive data collapse of the Binder cumulant as a function 
of Lr/L^ cf. Eq. (|4|. 

Before moving on, we comment on the two interre- 
lated (subleading) finite-size effects in the crossing point 
of Fig. [2} the crossing point between two subsequent 
Binder curves moves towards lower coupling for increas- 
ing system size, and accordingly the Binder cumulant 
at the crossing point decreases for increasing L. Con- 
sequently, the value of g* {K = Kc) will never be inde- 
pendent of system size L for finite systems. However, in 
our experience this vertical deviation from collapse of the 
Binder curves - which is particularly evident when focus- 
ing on the peak of the Binder curves as in our analysis - 
does not itself affect the finite-size estimate for z. More 
important is a possible horizontal deviation. Likewise, 
a slow convergence of the crossing points to Kc compli- 
cates the determination of the critical coupling for finite 
systems. The resulting uncertainty in z is dominated by 
this uncertainty in Kc, at least for the d = 2 case. 

It might be possible to obtain better precision for the 
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FIG. 2: (Color online) Illustration of the procedure for esti- 
mating the dynamical critical exponent z, as described in the 
text, here for the 2D quantum Ising model with site dissipa- 
tion and a — 0.2. a) The Binder cumulant y as a function 
of temporal system size Lt for a number of spatial system 
sizes L &t K = 0.160312. b) The peak value Binder cumu- 
lant g* as a function of coupling K. c) Finite-size analysis of 
the peak position of L* as a function of spatial system size 
L at criticality, Kc = 0.160312(2), which yields the estimate 
z = 1.97(3). The straight line represents a least squares fit to 
the data points. 



retained system sizes such that the value of z seems to 
have converged. For the case a = 0.2 considered above, 
the resulting estimate \s z — 1.97(3). No significant vari- 
ation in the dynamical critical exponent is observed for 
stronger dissipation, and we conclude that we have z ~ 2 
along the critical Une. However, we have not been able 
to determine conclusively whether or not one has ex- 
actly anomalous scaling dimension 77 = in the relation 
z = 2 — 77, which might be expectecP^ since the value 
d + z lies at the upper critical dimension for this phase 
transition for d — 2. 

We also give an estimate of the correlation length expo- 
nent v using the peak values g* (K) of the Binder cumu- 
lant. The leading order scaling properties of the Binder 
ratio can be stated afpH Q{K,L) = g{[K - ifji^/"), 
and assuming negligible finite-size effects in the obtained 
dimensions L*{L), one finds the finite-size relation 



log 



dK 



C 



1 



log L, 



(5) 



The slope dg*/dK is estimated by the finite difference 
Ag* over a small coupling interval around Kc, and C is an 
unimportant constant. The resulting finite-size analysis 
for a = 0.2 is illustrated in Fig. [3] and we find ly — 
0.49(1). This is very close to the expected (mean-field) 
value ly = 1/2 (Ref. fTU]). which is reasonable given that 
z « 2. 
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FIG. 3: Finite-size analysis for obtaining l/v for the 2D quan- 
tum Ising model. Here we have evaluated the slope of the 
Binder cumulant g* around K = 0.160312 for a = 0.2, which 
yields 1/ — 0.49(1). The straight line represents a least squares 
fit to the data points. 



critical coupling by using the finite-size analysis tech- 
nique presented in Ref. [30j for the crossing points, but 
in the present case with an additional (and unknown) 
finite-size effect in z, this more rigorous approach seems 
by no means straightforward. To ensure that finite-size 
effects are negligible, we have checked the dependence of 
z on the lowest value of L included in the fitting proce- 
dure. In the analysis illustrated in Fig. [2] we have only 



We finally note that, whereas increasing a does not 
lead to a significant change of z, it certainly does in- 
crease the prefactor a of the scaling relation Lr ~ 
and thereby the peak position L*. This reflects the in- 
creased anisotropy of the interactions, and can be seen 
also for a = when K and Kr are allowed to vary freely. 
At criticality one has a = 1 for Kt = K, with increas- 
ing a for increasing anisotropy K^/K. In fact, for the 
analytically solvable 2D Ising model there even exists an 
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exact mapping between system size anisotropy {i.e. 
and interaction anisotropy {i.e. Kr/K)W^ 
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III. QUANTUM ISING CHAIN WITH BOND 
DISSIPATION 

In this section, we will consider a (l+l)-diniensional 
quantum Ising model where the dissipative quantities of 
interest are bond variables involving Ising spins, rather 
than individual Ising spins themselves. The specific 
form of this dissipation kernel has been proposed as 
a candidate for describing the origin of the anoma- 
lous normal state properties of the cuprate high-T^ 
superconductors,'2Sl, but in that case involving two sets 
of Ising spin on each lattice point. Such a model, unlike 
the one we will consider, may be mapped onto a 4-state 
clock model, and may be approximated by an XY model 
with a four-fold symmetry breaking field, which in the 
classical case in two spatial dimensions is perturbatively 
irrelevant near criticality on the disordered side. Due to 
the degrees of freedom in our model being Ising spins 
with a spin gap, the present model should therefore not 
be regarded as directly comparable to a dissipative XY 
model that the authors of Ref. [5D] consider. It should 
rather be regarded as a simple, but spatially extended 
model system, illustrating how bond dissipation can af- 
fect a quantum phase transition, which is certainly an 
important question on its own right. 

In Fourier space the action is given by 



(6) 



The real space representation of this system is given by 
the action 



L 



S = —K ^ ^ <Jx,T<7x+l,T + Kr ^ ^ Ux,T'yx,T+l (7) 

X—1 T — 1 X — 1 T 



a 



cf. the site dissipation case in Eq. ([2]). Here, t^a^^r = 

The interpretation of this representation remains 
mostly the same as in the previous section. The only 
difference is that the coupling to the heat bath is given 
in terms of the Ising field gradients rather than the Ising 
fields themselves. In the limit q — 0, a; — > we may 
anticipate from the Fourier representation of the action 
that the last term becomes irrelevant, which implies the 
value z — \ for the dynamical critical exponent. It is also 
evident from Eq. Q that the bond dissipation is less 
effective than site dissipation in reducing quantum fluc- 
tuations. While site dissipation tends to align all spins 
in the temporal direction, the bond dissipation tends to 
align the difference of nearest-neighbor spins along the 



Trotter slices. At least in the presence of a finite coupling 
if 7^ this is a less effective way of reducing temporal 
fluctuations of individual spins than onsite dissipation. 

When expanding the dissipative term, it becomes clear 
that it contributes to both ferromagnetic and antiferro- 
magnetic long-range interaction. This renders the system 
intractable to the Luijten-Blote variant cluster algorithm 
used in the previous section. This algorithm builds up 
clusters with sizes comparable to the entire system and 
flips these as a whole, resulting in extreme correlations.'^ 
No cluster algorithm that effectively handles competing 
interactions has come to the authors' attention. 

In the Monte Carlo simu lations, we have therefore used 
a parallel temperin^^lM] algorithm which adequately 
handles the critical slowing down in the critical regime. 
A number of independent systems perform random walks 
in the space of coupling values, and this enables the sys- 
tems to effectively explore a rugged energy landscape like 
the one generated by the dissipation term in Eq. ([t]). 

The coupling values are distributed according to the 
iteration procedure introduced by Hukushima,ES! which 
renders the accept ratio of the attempted exchange of 
two adjacent coupling values independent of the coupling 
value. Consequently, the systems are allowed to wander 
relatively freely through the space of coupling values, al- 
though even more sophisticated distribution algorithms 
are available in that respect.!^ 




FIG. 4: Phase diagram of the quantum Ising chain with bond 
dissipation for AT^ — In (1 -|- \/2)/2. The ordered phase is 
found for large values of spatial coupling K and dissipation 
strength a. The filled square on the a axis represents an upper 
bracket for critical coupling Oc when the spatial coupling is 
tuned to zero, see the text. 

The parameter is fixed at ln(l + \/2)/2 w 0.4407, 
the critical coupling Kf^ is thus the same as for the 
isotropic 2D Ising model when the dissipation strength 
is tuned to zero. Anticipating z — \^ this choice also 
ensures that the simulations will be performed for con- 
venient values of L and i^. The further steps necessary 
to find information about the critical properties are the 
same as discussed in section |ll| The phase diagram of the 
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FIG. 5: (Color online) Data collapse of the Binder cumulant 
for z = 1 at Kc = 0.3306 for the ID quantum Ising chain with 
bond dissipation and a — 0.1. The error bars are obtained 
from a jackknife analysis in the reweighting procedure. Inset: 
Finite size analysis resulting in dynamical critical exponent 
z = 1.007(15). The straight line represents a least squares fit 
to the data points. 

system in the a-K plane is shown in Fig. [4] 

For this model, the critical exponents are extracted 
for the two dissipation strengths a = 0.1 and 0.2. In 
Fig. [5] we show the results for the dynamical critical 
exponent for a = 0.1 as illustrated by the collapse of the 
Binder cumulant curves discussed in Sec. Ullfor the value 
z = 1. The results confirm the proposed value of z based 
in naive scaling arguments, and it appears that the bond 
dissipation term is indeed irrelevant. 

The value of the dynamical critical exponent is very 
sensitive to finite size effects and therefore challenging to 
obtain with the algorithm we have used given the limi- 
tations this entails. Increasing the dissipation strength 
makes these challenges more apparent, so to illustrate 
the dependence of z on system size we plot in Fig. [6] 
z as a function of system size for a fixed K — Kc for 
a = 0.2. Note that three adjacent L values have been 
used to calculate every value for z, {L) denoting the av- 
erage of these. The evolution of z is clearly seen to ap- 
proach z 1 in the thermodynamic limit. Even larger 
dissipation strengths tend to require much larger system 
sizes not practically feasible with the current algorithm. 
Results for such dissipation strengths are therefore not 
included here. 

We have also attempted to extract the correlation 
length exponent v for both dissipation strengths. When 
discarding the smallest system sizes where finite size ef- 
fects are expected to be important, the values are found 
to be = 1.00(2) for a = 0.1 and f = 1.005(8) for 
a = 0.2. This corresponds well with the exact value 
1^ = 1 expected for the universality class of a 2D Ising 
model. 

Sufficiently strong dissipation brings the critical cou- 
pling Kc towards zero, and, as indicated on the a-axis 



FIG. 6: The evolution of z as a function of system size for the 
ID quantum Ising chain with bond dissipation with a = 0.2. 
Every point is calculated for the same coupling value Kc = 
0.231 from three adjacent system sizes, and the error bars are 
obtained from the least squares fit in the finite size analysis. 



of the phase diagram in Fig. [4| the model undergoes a 
purely dissipative phase transition at some critical dissi- 
pation strength ac- The ground state at K = consists 
of columns in the direction of imaginary time of ordered 
Ising spins. However, the direction of ordering is in gen- 
eral not uniform, as can be seen from Eq. ([7|), since a 
column can be flipped as a whole with no cost of energy. 
This nonuniform order prohibits the use of Binder cumu- 
lant curves to determine the critical coupling, so the ex- 
act value of ac is difficult to deduct from the simulations. 
These obstacles make an estimate of the dynamical crit- 
ical exponent unfeasible by our methods. Furthermore, 
since this phase transition is not of Kosterlitz-Thouless 
nature, any variety of the method of Ref. [17]also seem to 
be inapplicable to this model. 

To corroborate that there is in fact a phase transi- 
tion to an ordered state for increasing a also aX K = 0, 
we present in Fig. [7] results for the temporal spin-spin- 
correlation g^r) = {(Jx,T'^xfi) ■ It is clear that this corre- 
lation function decays exponentially to zero for low dissi- 
pation strengths, while in the opposite limit of strong dis- 
sipation the correlation function quickly decays to some 
finite value. The character of the correlation function as 
a is tuned through the intermediate region is better illus- 
trated in Fig. [8) where we have extracted the temporal 
correlation length ^7-. The diverging correlation length 
signifies a critical region with algebraic decay of the cor- 
relation function. The spatial correlation length ^, on 
the other hand, we have found to be vanishing also in the 
critical region, and the behavior of the system depends 
only very weakly on its spatial extent L. From a crude 
finite-size analysis based on Fig. [Sj we obtain the value 
ac < 0.64 as a best estimate for a upper bracket of the 
critical coupling, as we indicated in the phase diagram in 
Fig.|4} 
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FIG. 7: The temporal spin-spin correlation function (;(r) = 
{o'x,tO'x,o) for the ID quantum Ising chain with bond dissi- 
pation at = with L = 20 and Lt = 600. The decay of 
the correlation function is illustrated for four different dissi- 
pation strengths as the system goes from the disordered phase 
(a = 0.52) to the ordered phase (a — 0.68). 
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FIG. 8: (Color online) The temporal correlation length 
as a function of dissipation strength a for the ID quantum 
Ising chain with bond dissipation a.t K — 0. Because of the 
extremely anisotropic scaling in this limit, and to be sure 
to avoid any spatial finite size effects, we have chosen to fix 
L = 20. 



IV. DISCUSSION 

We will begin the discussion of our results by taking a 
closer look at the scaling arguments presented in Sec. |l] 
for finding the dynamical critical exponent. As indicated 
here, one important caveat of such arguments is that they 
only tell what exponent is naively expected to the low- 
est order approximation, and in general field-theoretical 
methods (see, e.g., Ref. fTO|) are needed to ascertain how 
higher order corrections modify this estimate. Further- 
more, with several terms in the quadratic part of the 



action, it is not always obvious which terms should be 
required to balance at the critical point, or for which 
phase transitions this is valid. 

For site dissipation, one obtains z = 2 by balancing 
the spatial term and the dissipative term, since the dy- 
namic term will be subdominant to the dissipative term 
for all positive z. For the bond dissipation case a similar 
argument excludes the possibility z — 2 for which the 
dissipative term would balance the dynamic term, since 
they both would be subdominant to the spatial term for 
all z > 1. It is therefore interesting to ask if the pos- 
sibility z = 0, or alternatively z ^ 1, can be consid- 
ered. In the limit that z is strictly zero, a dissipative 
term on the form |a;| would balance the dynamic term 
whereas a dissipative term on the form q^|a;| would bal- 
ance the spatial term, but in the latter case both would 
be subdominant to the dynamic term. One interpreta- 
tion is that 2 = in both cases would imply unrestrained 
quantum fluctuations resulting in spatial correlations be- 
ing infinitely stronger than temporal correlations, so that 
each Trotter slice is essentially independent. In this in- 
terpretation, a strictly vanishing dynamical critical ex- 
ponent may however be considered unphysical since we 
are assuming a transition to uniform order for the entire 
((i-l-l)-dimensional system by taking the limit q — >■ 0, 
a; 0. 

Likewise, tuning Kr — may be considered unphysical 
since one removes the origin of the quantum nature of the 
system. For this reason one can not say that there will 
exist a quantum phase transition with z = for the bond 
dissipation model even if the term had been removed 
from the action. The origin of the term in a physical 
quantum model can be a transverse magnetic field in the 
Ising case or a Josephson charging energy in the XY case, 
and the interpretation of the prefactor Kr is in general 
the inertia of the degrees of freedom. Even though we 
have chosen to operate with a nonspecific parameter K^-, 
we therefore do not regard taking = admissible in 
our simulations. 

The opposite limit of z = oo may similarly be inter- 
preted as spatially local criticality with correlations in the 
imaginary time direction independent of (the vanishing) 
correlations in the spatial directions, see, e.g., Refs. HOI 
[la and [m This is trivially the case in the limit K = 
for site dissipation with a > ac, although one may argue 
that z is undefined in that case as the system is strictly 
decoupled in the spatial directions. The same argument 
can not be applied to bond dissipation. For that model, 
the system does not experience dimensional reduction as 
K —i' 0, but is still dependent (although very weakly) 
on the spatial extent of the (d+l)-dimensional system. 
We should however note that the approach taken here 
for determining the dynamical critical exponent is not 
applicable when z is either strictly zero or infinite, and 
also for a constant value z 3> 1 it would be very difficult 
to determine the dynamical critical exponent for practi- 
cally attainable lattice sizes. If, on the other hand, one 
has z — oo in the sense of activated dynamical scaling, 
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the method is in principle feasiblel^Sl 

Before continuing the discussion of the bond dissipa- 
tion, we comment further on the relation between the 
real space representation of q^|a;| and the form of the 
bond dissipation used in Eq. Q. When Fourier trans- 
forming q^|w|crq.tjcr_q^_tj from Eq. (|6| and discretizing 
the resulting differential operators, we arrive at 

\Lr J sm [tt/Lt-\t - r'l) 

Now, writing out the last term of Eq. ([7| and comparing 
with Eq. Q shows that the Fourier space representation 
of the bond dissipation can be written as 

S'bond = (q^kl -I- C"q2)crq,„cr_q,_„. (9) 

Here, C" depends weakly on dimensions for finite sys- 
tems. In other words, the bond dissipation is of the same 
form as o^\lo\ dissipation, but with renormalized spatial 
nearest neighbor coupling, which however does not al- 
ter the critical exponents of the model. This extra term 
originates with the counterterm introduced to cancel out 
the renormalization of the bare potential that arises due 
to the coupling with a heat bath.l^l For the Ising model, 
this renormalization is responsible for stabilizing ferro- 
magnetic order at if > 0. 

We will now turn to the analysis of simulations on finite 
lattices, in particular with respect to the scaling relation 
Lt = aL^ and the system anisotropy expressed by it. 
To interpret our results it is useful to consider the de- 
pendence of both z and a on the dissipation strength a, 
and the variation of these quantities can be understood 
as follows. If the dissipation term is relevant and thus 
determining the universality class, we may assume that 
the value of z will be given by the form of this term even 
for infinitesimal a > in the thermodynamic limit. In 
this case, increasing the dissipation strength a further 
will therefore not change z, but the prefactor a will have 
to change to reflect the increased interaction anisotropy. 
Correspondingly, when the dissipation term is an irrele- 
vant perturbation, the dynamical critical exponent will 
remain z = 1 in the thermodynamic limit. Upon increas- 
ing a, the dissipation will never grow strong enough to 
alter the universality class, but the non-universal prefac- 
tor a will in general change also in this case, and whether 
it increases or decreases is determined by how the dissi- 
pation changes the overall interaction anisotropy. 

Regarding the evolution of a upon increasing a for 
the bond dissipation case, there are now two effects that 
must be considered separately. One implicit effect is that 
increasing a decreases K — at criticality, thereby 
increasing the anisotropy ratio K^/K, which results in 
a much larger a for large values of a. The other ef- 
fect is that arising explicitly from the dissipation term 
and its contribution to the effective coupling strength 
in the imaginary time direction. Whereas a site dissi- 
pation term obviously increases the anisotropy when in- 
creasing the dissipation strength while keeping the other 



coupling values fixed, such an enhancement of a does 
not appear for bond dissipation. This can be seen - as 
we have checked - by evaluating a for increasing a for 
isotropic short-range coupling, i.e. Kr = K. One pos- 
sible interpretation of this result is that although bond 
dissipation does not change universality, it favors z < 1 
behavior, which can also be recognized from Fig. [6j In 
other words, the dissipation term contributes to making 
the temporal dimension less ordered than the spatial di- 
mension, in strong contrast to the case of site dissipation. 
This would in part explain why one needs much longer 
simulations and larger systems to obtain reliable results 
for strong bond dissipation. 

Given that the exceedingly strong finite-size effects 
thwart a precise determination of z for higher values of a, 
one should in general also consider the possibility of con- 
tinuously varying critical exponents. However, we have 
shown that z w 1 for a = 0.1 and presented solid argu- 
ments favoring that this is the case also for a = 0.2, as 
it is obviously also in the limit a = 0. Therefore, if the 
exponents are in fact continuously varying, they begin to 
vary only for dissipation strengths above a > 0.2, and 
would furthermore have to be varying very slowly. 



V. CONCLUSION 

This work represents a further step towards simula- 
tions of physically interesting extended quantum systems 
with dissipation. Using Monte Carlo methods, we have 
studied a model similar to that by Werner et al. ^ but 
with higher spatial dimensionality, as well as a model 
with one spatial dimension but with bond dissipation in- 
stead of site dissipation. We have found that the (2-1-1)- 
dimensional model with site dissipation has a dynamical 
critical exponent very close to the corresponding d = 1 
model, i.e. z ss 2. Bond dissipation, on the other hand, is 
fundamentally different, and our results strongly suggest 
that this form of dissipation is irrelevant to the univer- 
sality class, i.e. z « 1 and non- varying. We therefore 
believe that the same dynamical critical exponent also 
applies to (2-|-l)-dimensional models with bond dissipa- 
tion for the same degrees of freedom, although we have 
not been able to reach sufficiently large systems to show 
this convincingly by numerical means. In both cases, the 
numerical estimates for the dynamical critical exponent 
is consistent with those found by naive scaling arguments 
on the quadratic part of the action. 
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